Accurate Weak Lensing of Standard Candles, Part 1: Flexible Cosmological Fits 
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With the availability of thousands of type la supernovae in the near future the magnitude scatter 
induced by lensing will become a major issue as it affects parameter estimation. Current N-body 
simulations are too time consuming to be integrated in the likelihood analyses used for estimating 
the cosmological parameters. In this paper we show that in the weak lensing regime a statistical 
numerical approximation produces accurate results orders of magnitude faster. We write down 
simple fits to the second, third and fourth central moments of the lensing magnification probability 
distribution as a function of redshift, of the power spectrum normalization and of the present-day 
matter density. We also improve upon existing models of lensing variance and show that a shifted 
lognormal distribution fits well the numerical one. These fits can be easily employed in cosmological 
likelihood analyses. Moreover, our theoretical predictions make it possible to invert the problem 
and begin using supernovae lensing to constrain the cosmological parameters. 
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I. INTRODUCTION 

More than a decade after the discovery of accelera- 
tion [1, 2], type la supernovae (SNe) are still the most ac- 
curate tool to map the cosmic expansion at low redshift. 
Intense search programs for the next years, like DES [3] 
and the LSST [4] will make available to the cosmologists 
thousands of new supernovae, enhancing the use of these 
standard candles as cosmology probes. With the honor 
of being the most useful cosmology probe comes the duty 
of understanding and reducing the impact of possible sys- 
tematic biases [5], for instance dust extinction, progenitor 
contamination, intrinsic variability etc. 

Among these systematic effect, perhaps the most in- 
teresting from a cosmological point of view is the magni- 
fication lensing. The light ray bundle from any source at 
large redshift in an inhomogeneous universe is of course 
subject to magnification or demagnification that mixes 
with the intrinsic source magnitude and change its statis- 
tical properties. The intrinsic scatter among supernovae 
at a given redshift is therefore convolved with the scat- 
ter induced by the matter fluctuations along the line of 
sight. Since the fluctuations themselves depend on the 
cosmological parameters, the lensing noise is actually an 
additional signal that can help with parameter estima- 
tion. We pursue this avenue in [6]. In the present paper 
instead we focus on estimating the lensing distortions 
through the second, third and fourth central moments of 
the lensing PDF. 

The lensing probability distribution function (PDF) 
is the fundamental quantity describing the statistical 
(de)magnification of distant sources by large-scale inho- 
mogeneities. After some pioneering work (e.g. Ref. [7]), 
various approaches have been followed so as to compute 
the lensing PDF relative to a given cosmology. A first 
approach (e.g. [8-11]) relates a "universal" form of the 
lensing PDF to the variance of the convergence, which is 
fixed by the amplitude ag of the power spectrum. More- 



over the coefficients of the proposed PDF may be trained 
on a grid of N-body simulations. A second approach 
(e.g. [12-18]) is to build a model of the universe, e.g. by 
means of an N-body simulation, and directly compute the 
relative lensing PDF, usually through time-consuming 
ray-tracing techniques. The flexibility of this method is 
therefore penalized by the increased computational time. 

In order to combine the flexibility in modeling with 
a fast performance in obtaining the lensing PDF the 
stochastic gravitational lensing (sGL) method has been 
introduced in Refs [19-21]. The sGL method is based on 
the weak lensing approximation and generating stochas- 
tic configurations of inhomogeneities along the line of 
sight. For instance, its speed allows one to include the 
full (cosmology-dependent) weak- lensing effects in the su- 
pernova analysis. This was first carried out in [22] for 
a simple toy-model, consisting of a universe populated 
by a distribution of large halos, all with the same mass, 
0(1O 14 M ). 

The outline of this paper is as follows. In Section II 
we will describe the universe model we adopt and briefly 
review the very basics of sGL. In Section III we will give 
our results as far as the lensing moments are concerned, 
while in Section IV we will show how to reconstruct the 
lensing distribution from the latter moments. Finally, 
we will conclude in Section V. In the Appendixes A and 
B, we discuss some properties of lensing and the errors 
intrinsic in our approach, respectively. 



II. SETUP 

A. Universe model 

We consider homogeneous and isotropic Fricdmann- 
Lemaitre-Robertson- Walker (FLRW) background solu- 
tions to Einstein's equations, on top of which we add 
matter perturbations describing virialized halos. We as- 



sume as fiducial a model sometimes referred to as "open" 
(non-flat) wCDM, for which the parameters describing 
the background are the present-day expansion rate h, 
matter density parameter £l m o, dark-energy density pa- 
rameter £1\q and dark-energy equation of state parame- 
ter w. 

We will model the matter contrast 5m (r, t) according 
to the so-called "halo model" (HM) (see, for example, 
[23-30]), where the inhomogeneous universe is approxi- 
mated as a collection of different types of halos whose 
positions satisfy the linear power spectrum. More pre- 
cisely, the halo model assumes that on small scales (large 
wave numbers k) the statistics of matter correlations are 
dominated by the internal halo density profiles, while on 
large scales the halos are assumed to cluster according to 
linear theory. In other words, the nonlinear evolution is 
assumed to produce only concentrated halos. The model 
does not include intermediate extended structures such 
as filaments and walls. See [20] for an improved halo 
model that also includes filaments. 

The parameters describing the spectrum of perturba- 
tions are, besides the ones relative to the background, the 
power spectrum normalization as, spectral index n s and 
baryon density parameter flbo- For the transfer function 
T(k) we use the fit provided by the Equations (28)-(31) 
of Ref. [31]. See [20] for more details about the power 
spectrum used. 

Regarding the modeling of halos, we use the halo mass 
function given in Eq. (B3) of Ref. [32], which is defined 
relative to a spherical-overdensity halo finder and has a 
good degree of universality [33]. The halo profiles are 
modeled according to the Navarro-Frenk- White (NFW) 
profile [34] , which is able to model both galaxy-sized ha- 
los and superclusters with a appropriately chosen concen- 
tration parameter c(M, z). The concentration parameter 
depends on the cosmology and we use the universal model 
proposed in Ref. [35] , which allows for an accurate deter- 
mination of a halo's concentration by relating the latter 
to the time at which the halo's main progenitor first as- 
sembled 4% of its final mass. We however neglect corre- 
lations in the halo positions. As shown in [20, 21], this 
should be a good approximation for the redshift range of 
z < 1 in which we are mainly interested in this paper. 

Summarizing, our universe model, which is determined 
by the vector of parameters {h, fl mQ ,fl\ , W,ag, n s , Qbo}, 
approximates in a consistent way the real universe at the 
nonlinear scales of clusters, the ones relevant for lensing 
of standard candles. 



B. Lensing model 

As we will now see, the sGL method predicts the sta- 
tistical distribution of the lens convergence k, which in 
the weak-lensing approximation is given by the following 
integral evaluated along the unperturbed light path [36] : 



where the quantity 5m{Ti t) is the local matter density 
contrast, the density puc = «o Pmo is the constant mat- 
ter density in a co-moving volume, and we defined the 
auxiliary function 



G(r,r s ) 



^G f k (r)f k (r s - r) 
c 2 a fk(r s ) 



(2) 



which gives the optical weight of a matter structure 
at the comoving radius r. The functions a(t) and 
t(r) are the scale factor and geodesic time for the 
background FLRW model; r s — r(z s ) is the comov- 
ing position of the source at redshift z s ; and finally 
fk{r) = sm(rvk)/vk, r, sinh(r%/— k)/\/— k depending 
on whether the curvature k is >, =, < 0, respectively. At 
linear level, the shift in the distance modulus caused by 
lensing is expressed in terms of the convergence only: 



Am(z) ~ 51og 10 [l - k(z)] ~ - 



log 10 



k{z) . (3) 



k(z s ) 



dr p AI c G(r, r s ) 5 M (r, t(r)) 



(1) 



Eqs. (1) and (3) illustrate how for a lower-than-FLRW 
column density the light is demagnificd, while in the op- 
posite case it is magnified. Also note that Eq. (1) is only 
valid for point sources (smoothing angle 8 — 0). 

Eq. (1) connects the statistical distribution of mat- 
ter to the statistical distribution of convergences. The 
sGL method for computing the lens convergence is based 
on generating random configurations of halos along the 
line of sight and computing the associated integral in 
Eq. (1) by binning into a number of independent lens 
planes. Because the halos are randomly placed their oc- 
cupation numbers in parameter-space volume cells follow 
the Poisson statistics. This allows rewriting Eq. (1) as 
a sum over these cells characterized by the correspond- 
ing Poisson occupation numbers. The various individ- 
ual contributions to the convergence are then additively 
combined. By generating many halo configurations one 
can sample the convergence PDF. A detailed explana- 
tion of the sGL method, which can also model filamen- 
tary structures confining halos, can be found in [19-21] 
and a publicly-available numerical implementation, the 
turboGL package, at turbogl.org. We note that although 
turboGL is not strictly speaking a ray-tracing procedure 
it emulates it by constructing many halo realizations 
along a given direction. To obtain the results in this 
paper we ran turboGL with a statistics of 1.5 • 10 6 "light 
rays". In other words, we shoot 1.5 million light rays for 
every point in the cosmological parameter space. 

Overall, we estimate that our results can be relied upon 
at the level of ~10%. We discuss in more detail in Ap- 
pendix B the full error budget regarding this estimate 
and the approximations of this setup. 



III. MOMENTS OF THE LENSING PDF 

The modeling described in the previous Section and 
implemented in (the latest) version 2.0 of turboGL, al- 
though many orders of magnitude faster than running 
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Figure 1. Comparison of the convergence PDF's from turboGL [19, 20] and ray-tracing from two N-body simulations, dubbed 
H07 [15, 16] and Til [17]. The former is based on the Millennium Simulation (at z = 1.08), the latter on a more recent 
simulation (at z = 1). Left: slight discrepancy at negative k; although the cosmological parameters are different on both 
simulations, the discrepancy is relevant. Right: zoom on the high-magnification tail and comparison of image and source plane 
PDF's. See Section III A and Appendix A for more details. 



an N-body simulation and relative ray-tracing routine, 
still needs about 15-20 seconds to generate a PDF with 
a statistics of 10 5 using one core of a CPU at 2-3 GHz. 
Therefore, it can be computationally expensive to run 
sGL on a grid or Monte-Carlo for many points. It is 
therefore convenient to create shortcuts which could al- 
low for a fast estimation of the weak lensing distortions 
in different cosmologies. In fact, as discussed before, dif- 
ferent shortcuts to the full problem were proposed in the 
literature in recent years, see e.g. [8-11]. These, how- 
ever, usually take the form of a simplified algorithm to 
compute the lensing PDF, the application of which is not 
completely straightforward and still adds some computa- 
tional complexity to any likelihood estimation. 

Here we take a different approach. Instead of giving 
a prescription so as to obtain the lensing PDF, we focus 
on the information condensed into its central moments 
M2-4) which are defined as 



tx n = ((x-(x)r), 



(4) 



and are directly related to the variance, skewness and 
kurtosis of the PDF. We provide below analytical fitting 
functions for these three moments with respect to the 
most relevant cosmological parameters, to wit er 8 and 
fi m o, as we U as redshift, all for the very broad ranges 
listed in (5). 

There are a number of advantages in this approach. 
First, as will be shown, it makes it very straightforward 
to include lensing in any likelihood analysis. Second, it 
adds almost no computational complexity. Third, the 
moments may be directly related to observables. For in- 
stance, if the supernovae are intrinsically gaussian, the 
skewness of the lensing PDF would be exactly the skew- 
ness of the final convolved distribution. In other words, 
our estimation of the skewness at a given redshift would 
be the skewness one would measure in the distributions 



of standard candles in a thin bin at that redshift. This 
in turn also allows one to quickly get answers to ques- 
tions regarding what might be called the inverse prob- 
lem of standard candle lensing: how can a measurement 
of the actual supernova PDF (or of any other standard 
candle candidate) in different redshift bins constrain cos- 
mological parameters such as as and fl m o- We analyse 
this latter prospect in detail in a forthcoming companion 
publication [6], and instead focus here on the moments 
themselves and how to use them. 



A. Comparison with previous results 

We will now compare the PDF generated by turboGL 
with previous results obtained by shooting rays through 
the matter field of N-body simulations. In particular we 
will consider the findings of Refs. [15, 16] (dubbed H07) 
which are based on the ACDM model of the Millennium 
Simulation [37] consistent with the 1-year WMAP re- 
sults [38] and the ones of Ref. [17] (dubbed Til) which is 
based on an N-body simulation consistent with the more 
recent 5-year WMAP results [39]. 

Figure 1 shows an overall agreement between the con- 
vergence PDF's obtained with turboGL and through ray 
tracing. The left panel shows a very good agreement be- 
tween the minimum magnification of turboGL and H07. 
The agreement is somewhat worse for the case of Til, 
which gives a minimum magnification lower than both 
turboGL and H07. Regarding H07, this discrepancy 
could be due to the different smoothing adopted of the 
matter field [40]. Regarding turboGL, it could be due to 
the use of the halo model for the matter field; indeed, as 
shown in [20] , the minimum magnification shifts to more 
negative values if one includes filamentary structures. It 
thus seems that for very small convergences, the differ- 



ent limitations of turboGL and H07 produce the same 
results. 

We should stress at this point that while turboGL is 
limited to the weak-lensing regime, the analyses of H07 
and Til consider full nonlinear lensing. As discussed 
in the Appendix A, one can compute the distribution of 
lensing distortions in either the source (S) or image (I) 
plane. The reference plane used by the sGL method and 
thus by turboGL is the S plane, and we will focus on 
that plane since it is the one in which lensing is directly 
related to observed quantities. Strictly speaking, in the 
deep weak lensing limit of k <C 1 this distinction is irrel- 
evant. However, as we will consider convergences up to 
K cut — 0.35, the output of turboGL should be compared 
to results relative to the S plane. The difference between 
the PDF's in the two planes is substantial only in the 
high-magnification tail, as shown in the right panel of 
Figure 1. There it is shown how turboGL agrees remark- 
ably well with both Til and H07 up to a K cu t = 0.35. 
In fact, the noisy oscillations at high convergences can 
be made better by hust using higher statistics (here we 
plot results with the statistics employed throughout this 
work - see Section II B). If one considers the original HO 7 
result (which was computed in the image plane) a small 
disagreement between turboGL is found. In summary, 
for small to medium convergences, turboGL is in very 
good agreement with both Til and H07, as long as one 
considers the source plane. 

In this paper we are mainly interested in the moments 
of the lensing PDF. We will now test, using Til as bench- 
mark, the performance of turboGL as far as the moments 
are concerned. We will not further consider H07 as we 
will restrict our analysis to the lensing PDF in the source 
plane. Figure 2 shows the comparison of turboGL and 
Til for the second-to-fourth central moments ^2-4,icns 
at the two source redshifts of z = 1 and z = 2. The plots 
are with respect to K C ut, the convergence above which 
we set the PDF to zero. The plots show a very good 
agreement for z = 1. The agreement gets slightly worse 
at z = 2: this could be due to the fact that the halo 
model works less well at the higher redshifts at which 
non-virialized structures (neglected in the halo model) 
play a more important role. However, since most super- 
novae lie in the z < 1 region, we can rely on turboGL 
to very good accuracy. Based on these results, we will 
adopt in the following the value K cu t = 0.35. 

This can be taken as a somewhat conservative value, 
since based on Figure 2 higher values seem to be in very 
good agreement for z < 1. We nevertheless want to pro- 
vide estimates for higher redshifts and, more importantly, 
for different cosmologies, for which N-body ray tracing 
lensing PDF's are currently not available. Although we 
do not expect a substantially different agreement for dif- 
ferent values of f2 m o or <t 8 , we currently cannot test this 
with N-body simulations. 

Finally, we note that one in principle could try and im- 
plement a K cut on the real data, by excluding the (very 
rare) standard candles for which there is independent evi- 
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Figure 2. Comparison between the N-body lensing PDF 
of Til [17] (in the source plane) and turboGL for different 
choices of K cut and for both z = 1 and 2 = 2. K cut is the 
convergence above which we set the PDF to zero. The agree- 
ment is very good for z = 1 but gets slightly worse for the 
higher moments at higher redshifts, especially for /13. Since 
most supernovae lie in the z < 1 region, we conclude that we 
can rely on turboGL to very good accuracy. We have chosen 
K cu t = 0.35 for our fits. All plots are in units of magnitudes, 
see Eq. (3). See Section III A for more details. 



dence of either strong or "medium" lensing. In fact, using 
the full lensing PDF of Til we estimate that less than 
1 in 20000 supernovae at z = 1 should have K cu t > 0.35 
(and thus magnification larger than 1.8), so a more care- 
ful check of high magnification candidates might be feasi- 
ble (see for instance [42]). Incidently, there was a recent 
claim in the literature [43] of a supernova with a huge 
magnification ^30; the probability for this to happen 
due to halos is less than one in a million, and in fact 
the authors suggest that the lens might be a compact 
object such as a black hole. An exclusion of high mag- 
nification events is thus important if one wants to use 
the higher moments, as they are the most sensitive to 
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Figure 3. Second-to-fourth central moments /ia-4, ions of the convergence PDF as a function of redshift and five different 
cosmological parameters {ag, Q m o, f^Ao, w , n s}, changing one at a time while keeping all the others fixed at the 9-year WMAP- 
only flat-ACDM best-fit values [41] and z = 1. As can be seen, apart from z it is clear that <jg and f2 m o are the most relevant 
quantities, whilst the dependence on ^ao, w an d n s is almost negligible. A value of K cu t = 0.35 has been adopted, and all plots 
are in units of magnitudes, see Eq. (3). See Section IIIB for more details. 



these. We nevertheless do not investigate further these 
complications in this work. 



B. Dependence on cosmological parameters 

In this Section we will quantify how strongly lens- 
ing depends on the relevant cosmological parameters. 
Figure 3 shows the change in the second-to-fourth cen- 
tral moments ^2-4,icns as a function of redshift and five 
different cosmological parameters {as,^l m o^Ao,w,n s }, 
changing one at a time while keeping all the others 
fixed at the 9-year WMAP-only fiducial values [41] and 
2 = 1. As we described above, we are assuming a non-flat 
wCDM model throughout this paper. It is clear from Fig- 
ure 3 that eg and f2 m o are the dominant parameters, con- 
sistently producing variations ~10 larger in observation- 
ally comparable ranges. This is what one would expect, 
as lensing is mostly governed by the amount of matter 



and related clustering. The dependence on £Iao, w and 
n s (and on h and Qbo, which were also computed but 
are not shown in Figure 3) is almost negligible. Based 
on these findings we decided to consider in the following 
only the dependence of lensing on z, as and fi m rj. 

C. Fitting functions for the moments 

Here we give analytical fitting functions for the second- 
to-fourth central moments /X2— 4,iens a s a function of 
{z, eg, ^mo} which are very accurate in the domain 



< z < 3, 

0.35 < er 8 < 1.25, 

0.1 < flmo < 0.52. 



(5) 



In fact, in this entire domain, the average RMS error is 
only 4% for all three moments. Using magnitudes (see 
Eq. (3)), the fitting formulae are: 
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Figure 4. Dependence of the weak-lensing standard deviation (in magnitudes, see Eq. (3)) with respect to {z, as, Q,mo} from 
the numerical results (red dots) and the analytical fitting functions (solid surfaces/lines). Top: 3D plots covering the full range 
of validity of the fitting functions. Bottom: same but depicted in 2D plots, with the addition of the fits found in [14] (gray 
dashed line) and [44] (gray dotted line). Note that our fits approximate the numerical results very well in the whole parameter 
space and are in good agreement with [44]. 



Clens(£,0- 8 ,Q m o) — 

l / 3 I n ^ _ 

^4 1cns( Z ' Cr 8i^mo) = 



0.0004 - 0.00176cr 8 + (-0.035 + a 8 n m0 + 0.0453a 8 )z _ 
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0.018+O m0 cT8Z 



0.3z 



(0) 
(7) 

(8) 



Figure 4 depicts the fits for ai cns for different parame- 
ters. The first row cover the full range of validity of the 
fitting functions. The middle row represents the same 
but depicted in 2D plots, with the addition of the Holz 
& Linder fit of am, — 0.093 z [14] (gray dashed line) and 
of best fit to the real supernovae data ajio — 0.055 z 
found by Jonsson et al. [44] (gray dotted line). Figure 5 
similarly shows the fits for /j, 3 and /i 4 . As can be seen, 
all three central moments are of the same order of mag- 
nitude. For all three moments the fit are almost always 
within 5% of the numerical values and have an average 
RMS error of 4% (see Appendix B for more details) . 

We would like to point out that the convergence PDF 
has also a small but nonzero mean. Although by con- 
struction the sGL method produces convergence distri- 



butions with zero mean (see Eq. (1)), the latter can be 
approximated (in the source plane) by 



(«) ~ -2( K 2 ) ~ - 
see Appendix A for more details. 



-2a 



lens ) 



(9) 



D. On the lensing variance 

We turn our attention now to our estimate of the lens- 
ing variance, Eq. (6). We find quite smaller values as 
compared to the cosmology-independent Holz & Linder 
fit of (Thl = 0.093 z [14], which although outdated is still 
used for example by the Supernova Cosmology Project 
(responsible for the Union Compilations - see [45]). This 
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Figure 5. Same as the bottom plot of Figure 4 for the case of the third and fourth moments. Once again, the fitting functions 
approximate the numerical results very well in the whole parameter space. Note that the use of fj, n instead of just /i n makes 
it clear that the first moments all have the same order of magnitude. 



is explicitly shown in the second row of Figure 4. Our 
results, besides being in very good agreement with Til, 
are also in better agreement with the observational con- 
straints obtained with real supernovae [44] (although we 
point out that the error bars render their results statis- 
tically compatible with [14]). This has important conse- 
quences that we detail below. 

In [14] they estimate that due to the extra scatter of 
lensing supernovae quickly loses "cosmological constrain- 
ing power" at high redshifts. For instance, assuming an 
intrinsic dispersion <7i nt = 0.10 mag, 3 supernovae would 
be needed at Z = 1.5 to provide the equivalent informa- 
tion of 1 supernova if there was no lensing. Here, we find 
that this "degradation factor", defined as 



lens 
°tat 



(10) 



is much smaller. Assuming fi m o = 0.3 and CTg = 0.8, only 
1.5 supernovae are needed to get the same information 
than what would otherwise be the case without lensing. 
If supernovae have larger intrinsic dispersion, say ai n t — 
0.15 mag, than the degradation is as low as 1.2. Figure 6 
depicts this for the intermediate value of a mt = 0.12 mag, 
for the parameter values given by WMAP9 [41], together 
with the allowed region by varying the parameters by 2er. 
We stress that this is very significant, as supernova ob- 
servational strategies and survey-planning may be over- 
neglecting the power of high-z candidates. Moreover, this 
alleviates one of the major issues for high-z standard can- 
dle candidates, such as radio galaxies [46], GRBs [47, 48], 



AGNs [49] or HII galaxies [50], all of which can be ob- 
served at z > 3. Since our fits fully take into account 
the dependence of lensing on cosmology, it allows a more 
careful estimate of the possible constraints of high-z stan- 
dard candles. We encourage, therefore, the scientific 
community to use our new fits for their own analyses, 
especially the fit of the lensing variance. 



IV. RECONSTRUCTING THE PDF FROM THE 
MOMENTS 

Here we will show how to reconstruct the lensing PDF 
from the moments provided by Eqs. (6-8). We will make 
use of the known fact [51] that the convergence distri- 
bution is well described by the lognormal PDF. This 
stems from the fact that the one-point distribution func- 
tion of the matter density field is itself well described by 
the lognormal distribution [52, 53] and that the conver- 
gence k is directly connected to the matter field 5m as 
shown by Eq. (1). The approach we will now take differs 
from [11, 51, 54]. We will now explain it in detail. 



A. Lognormal distribution 

The lognormal distribution depends on the value of two 
parameters, ^ gau and <T gau , which are the mean and the 
dispersion, respectively, of the gaussian distribution from 
which the lognormal is derived. The lognormal PDF, 



o- int = 0.12 mag 




Figure 6. Standard candle lensing degradation factor as a 
function of redshift. Green: our estimate for the best- fit pa- 
rameters of WMAP9, together with the allowed region by 
changing the parameters by 2a. Dashed gray: previous es- 
timate found in [14]. Dotted gray: previous estimate found 
in [44]. Big brown dots: numerical N-body results found in 
Til [17]. An intrinsic scatter of o"i nt = 0.12 mag is assumed. 



together with its mean and variance, can be written as: 
exp 



Jlogn\ x ) 



Mlogn 



(loga-jjgauj 

2<t!„ 



'2tt a 



gau X 

/2 



' logn ° 



l 



(ii) 

(12) 
(13) 



where x,<7 gau > 0. One can invert Eqs. (12)— (13) so as 
to express ^ gau , c g au as a function of the actual mean and 
dispersion of the lognormal distribution: 



Mgau = ^ lo S 

^L = log (i 



Mlogn 



1 + "logn/Mlogn 
i/Mlogn) • 



2 

°logn/ 



(14) 
(15) 



The lognormal distribution of Eq. (11) has a nonzero pos- 
itive mean, while the lensing distribution has an almost 
negligible mean R. So as to model the convergence PDF 
we will then translate Eq. (11) into our final template 
PDF for the convergence distribution: 



fn(x) 



exp 



[log(a:+/ii 



~ K )~ MgauJ 



2<tL 



27T Cr gau (x + Ailogn ~ «) 



(16) 



where ^ ga u, o" gau are given by Eqs. (14)-(15) as a function 

01 /-Mogn j ^logn ■ 



B. Using information from the moments 

The convergence PDF given in Eq. (16) depends on 
three parameters, the mean R, the dispersion <7i ogn and 
the shape parameter fi\ ogn (as the mean has been set, 
Miogn effectively constrains the higher moments of the f K ). 
The dispersion oi ogn can be simply fixed using Eq. (6): 



^logn — ^lcns • 

The mean can either be set to 

K=0, 

or to the (small) value given by 



n 2 
" z "lens ' 



(17) 



(18) 



(19) 



so as to take in consideration some higher order correc- 
tions, see Eq. (A5) and related discussion. 

We are then left with the shape parameter /ii ogn , which 
can be used either (i) to fix the third moment of f K to 
the ^3 lens of Eq. (7) or (ii) to fix the fourth moment of 
f K to the ^4 ! i ns of Eq. (8). We need then to express 
the shape parameter /ii ogn as a function of the third and 
fourth moments of f K . This can be done for case (i) and 
case (ii), respectively, by solving the following equations 
with respect to /Ji osn '- 



r 4 (3 + r 2 ) 



M3,lcns 

1^1 = r 4 (3 + 16r 2 + 15r 4 + g r 6 + r 8) 
Mfogn 



(20) 

(21) 



where r = <7i ogn //ii ogn . The previous equations can 
be easily solved analytically. The algebraic expression 
of the solutions, though straightforward to obtain with 
e.g. Mathematica or Maple, are however rather long and 
not particular illuminating. Consequently, they will not 
be explicitly reported here. Alternatively, we will give 
ready-to-use fitting functions. The solutions to Eqs. 
(20)— (21) depend on <7i cns and M3-4,lensj which in turn 
depend (mainly) on {z, ag, fi TO o}- It, therefore, natural 
to use directly the latter parameters. Fits to the solutions 
of (20)-(21) with RMS errors better than 4% within the 
ranges of (5) are, respectively: 
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Figure 7. Comparison of the reconstructed PDF's against the turboGL output. Left: the overall agreement is good as far as 
the body of the PDF is concerned. Right: the zoom on the high-magnification tail also shows a good level of agreement. The 
numerical values of the parameters specifying the convergence PDF f K of Eq. (16) are given in Table I. See Section IV for more 
details. 



where a = log 10 z, b — log 10 as and c = log 10 f2 



;io "mO- 



C. Example 

As an example of the formalism developed in the pre- 
vious sections we will consider the lcnsing PDF relative 
to the fiducial cosmology given by the 9-year WMAP- 
only flat-ACDM best-fit values [41] and z = 1. In Fig- 
ure 7 we compare the reconstructed PDFs against the 
turboGL output. We considered three cases: (i) using for 
the shape parameter the value given by (20) (or equiva- 
lently (22)) so as to match the third moment, (ii) using 
the value given by (21) (or equivalently (23)) so as to 
match the fourth moment and (iii) using the intermedi- 
ate value given by the mean of the last two. This makes 
sense as the two values of /zi ogn are numerically close to 
each other. To quantify the performance of the recon- 
structed PDF as far as the moments are concerned, we 
list in Table I the error in the moments as compared 
to the turboGL output. It is clear that f K is perform- 
ing well, particularly so in case (iii) where the errors are 
below the 2% level, well within the modeling errors of 
~10% intrinsic in our approach (see Appendix B). Also, 
even though the algebraic implementation of the lognor- 
mal distribution presented here differs from e.g. [54], the 
overall shape of the obtained distribution agrees rather 
well. 



V. CONCLUSIONS 

Next generation supernova searches will collect thou- 
sands of SNIa lightcurves up to high redshift, vastly ex- 
tending the present catalogs. In order to fully exploit 
this data abundance, the so-far almost negligible error 



criogn Miogn Error in ^3,lens Error in /tM.iene 



Case 



(i) 


0.0208 


0.0225 


0% 


3.8% 


(ii) 


0.0208 


0.0237 


-2.3% 


0% 


(iii) 


0.0208 


0.0231 


-1.2% 


1.8% 



Table I. Numerical value of the parameters specifying the con- 
vergence distributions plotted in Figure 7. The value of the 
shape parameter /ii ogn for case (iii) has been obtained by av- 
eraging the values relative to case (i) and (ii), which were 
obtained by matching the third and fourth moments of f K to 
the lensing moments given by turboGL. 



induced by lensing will need to be taken into account to 
avoid biasing the results. On the other hand, the lensing 
effect itself will become a source of cosmological and as- 
trophysical information, being related to the distribution 
and density profile of line-of-sight halos. 

The crucial step in promoting lensing scatter from 
noise to signal is to predict reliably its statistical prop- 
erties, quantified through their magnification PDF. The 
calibration of the lensing PDF through TV-body simula- 
tions is accurate but extremely time consuming and can- 
not, at this time, be carried out for a sufficiently large 
set of cosmological models as the one needed in the re- 
construction of the SNIa likelihood. 

In this paper we use the fast computation allowed by 
the sGL method [19-21] - which distributes halo and 
sources and evaluates the resulting lensing effect in the 
linear regime in a matter of seconds, i.e. several orders 
of magnitude faster than any ./V-body approach - to re- 
construct the lensing PDF for an array of cosmological 
models. First, we have shown that our fast results match 
those from state-of-the-art iV-body simulations. Then, 
we illustrated that the PDF is practically insensitive to 
h, £Iao,w, n s , flbo, while strongly dependent on O8,fi m o, 
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in agreement with expectations. We have quantified the 
latter dependence by plotting the behavior of the second, 
third and fourth central moments. We then provided ac- 
curate fits that give the moments for any redshift and any 
<78,f2 m o within a reasonable range. These fits have the 
advantage of being directly employable in a parameter 
estimation likelihood analysis. We quantify the error of 
the fits with respect to the TV-body simulations to be at 
the ~10% level. We also find that a shifted log-normal 
PDF - of which we give an easy-to-use prescription - 
approximates relatively well the numerical PDF. 

Our analysis considerably improves upon the popular 
fit [14] . In particular, we find that our estimate of the ad- 
ditional variance induced by lensing is quite smaller than 
in Ref. [14]. The lensing noise is found to be comparable 
to the intrinsic scatter at redshift z sw 3 rather than at 
Z « 1.3 as in Ref. [14], for a WMAP9 cosmology. This 
not only is in better agreement with preliminary observa- 
tional constraints of the lensing PDF [44] but also makes 
high-redshift supernova observations more relevant for 
cosmology and may entice reviews of supernova obser- 
vational strategies. 

It is clear that the next step in this line of research 
is to invert the problem and ask how well one can use 
the SNIa PDF to estimate the cosmological parameters, 
in particular a & ,fl m0 . If one assumes the SNe magni- 
tudes to be intrinsically Gaussian, then the entire ob- 
served third order moments and the excess-over-Gaussian 
fourth order moment are a result of lensing. If the rare 
strongly non-linear lensing events can be identified and 
eliminated from the analysis then the observations can 
be directly compared to our results and employed to con- 
strain (78, flmOi directly from lensing. On the other hand, 
if we fix the cosmological model and the observations 
show a significant deviation from the predicted lensing 
effect, then one would have an evidence for an intrinsic 
non-Gaussianity of the sources or a violation of lensing 
linearity. This work will be carried out in a subsequent 
paper. 



proven to be an exquisitely good one in all cases of phys- 
ical interest for cosmology. In such approximation, it is 
imperative to establish which is the plane of reference: 
the source plane or the image plane. In fact, the statisti- 
cal properties of the distributions of weak-lensing distor- 
tions vary significantly depending on this choice, as we 
will discuss below. In a broad sense, the source plane is 
closely related to observable quantities, whereas the image 
plane is often more directly related to numerical simula- 
tions. We shall denote by a subindex "S" ("I") averages 
conducted on the source (image) plane. 

A classic and general result of gravitational lensing 
is that photon conservation implies that the magnifica- 
tion (/2) PDF has unitary mean. This is true even outside 
the weak lensing regime, but this result assumes implic- 
itly that one is working on the source plane. In other 
words, (/2)s = M (fi)i- 

Since magnification /2, convergence n and shear 7 are 
related by 



I' 



one has for small k 



1 

(1-k) 2 -7 2 



(Al) 



(j, ~ 1 + 2k + 3k 2 + 7 2 . 


(A2) 


nplies that 




(/2) s = 1~1 + 2( K ) s + 3(k 2 ) s + (7 2 )s 


(A3) 


.-. ( K ) s c4(K 2 ) s -i( 7 2 )s. 


(A4) 



Now, in the weak-lensing approximation one has (n 2 )s = 
(7 2 )s [17]. Since turboGL does not compute the shear, 
we made use of the PDF's derived in [17] to test this 
numerically and found that, at z = 1, (k 2 )s — 4.7 1CP 4 
and (7 2 )s = 5.3 10 -4 , an agreement to the 10% level. 
Using this equality, the above simplifies to 



(*)t 



-2(k 2 ) s . 



(A5) 



ACKNOWLEDGMENTS 

It is a pleasure to thank Stcphan Hilbert, Kimmo 
Kainulainen, Eric Linder, Martin Makler, Bruno Moraes, 
Ribamar Reis, Ignacy Sawicki, Brian Schmidt, Peter 
Schneider and Ryuichi Takahashi for fruitful discussions. 
LA and VM acknowledge support from DFG through the 
TRR33 program "The Dark Universe". MQ is grateful to 
Brazilian research agencies CNPq and FAPERJ for sup- 
port and to ITP, Universitat Heidelberg for hospitality 
during part of the development of this project. 



Appendix A: Properties of the weak-lensing PDF 

In this work we are implicitly making use of 
the multiple-lens-plane approximation, which has been 



This correlation between mean and variance was indeed 
observed explicitly in [17] (see their Figure 5). 

If one repeats the above calculations in the image 
plane, results differ significantly. In the image plane there 
is a bias towards higher magnification events. This im- 
plies that 

(k)j = (kh) s =* (k)s + 2(k 2 ) s + 3(« 3 ) s + (k 7 2 )s . (A6) 

Using (A5) and up to higher order corrections, one finds 
that 



Wi 



0((k 3 



< K K )s 



(A7) 



Using once again the PDFs derived in [17] we find that, 
at z = 1, (k)i ~ — 0.005(k 2 )i, which should be confronted 
with (A5). 

Another important aspect of the weak-lensing PDF is 
the theoretical prediction of how the PDF of jl behaves 
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O m0 = {0.15, 0.3, 0.45}; cr 8 = {0.6, 0.8, 1.0} 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 



Figure 8. Percentage differences between fits and numerical 
results, showing that the fit is almost always within 5% of the 
numerical values, in agreement with the quoted average RMS 
error of 4% in the whole domain of Eq. (5). {Full, Dashed, 
Long-dashed} lines stand for Q m o = {0.15, 0.3, 0.45}, while 
{Blue, Green, Orange} lines stand for as = {0.6, 0.8, 1.0}. 



for /j> 1. Catastrophe theory predicts that, for a point 
source, in the high magnification limit one has [15, 55] 



PDF(/i)s ex M~ 3 : 
PDF(/z)i oc /i~ 2 . 



(A8) 
(A9) 



This in turn implies that even in the source plane, for 
a point source all central moments above 1 would be in- 
finite. In other words, the extra variance in a standard 
candle due to lensing would be unbounded. For an ex- 
tended source this is no longer the case and the PDF is 
suppressed for p, > /i max - Now jU max is typically well be- 
yond the weak-lensing validity range, even for galactic- 
scale sources. It was shown in [17] that for a circular 
source of radius 10 kpc/h at z = 5, /2 max — 4; for a su- 



pernova /imax must be much higher. Thus if one wants 
to make use of the weak-lensing central moment predic- 
tions, one must enforce a lower cut /ijjj™ to the PDF, 
which can in practice be achieved by removing all sources 
for which there is strong evidence that ft > /ijjj™ (see 
Section III A). 



Appendix B: Error budget of our results 

Our results are subjected to modeling errors. We 
we will quantify them as far as ui ens (= Aliens)' Aliens 

and /i 4 ' lcns are concerned. The first source of error 
comes from the weak-lensing approximation, which we 
estimated through comparison with exact results to be 
of about 5%. The second source is the matter model 
adopted. The halo model neglects indeed extended struc- 
tures such as filaments and walls. Based on the results 
of [20], we estimated this error at the 5% level. There 
are then the errors inherent to the mass function [32] and 
concentration parameter model [35] used. By changing 
the latter within the error ranges stated in the respective 
papers we found that they add a 3% error each to the 
overall error budget. Finally, there is the 4% error due 
to the fitting functions. Figure 8 gives the percentage 
differences between fits and numerical results, i.e. 



A|#» (% ) 



1/n, fit 



t/n, numerical 



(l 



V 



n. numer 



;ical 



(Bl) 



showing that for all three moments the fits are almost 
always within 5% of the numerical values. 

We list in Table II a summary of the errors intrinsic in 
the present modeling, which together do not exceed the 
10% level. These figures mainly refer to the second and 
third moment; the fourth nroment is subjected to a larger 
error of the order of 10 — 15%. Finally, though the fitting 
functions are accurate for the full ranges given in (5), 
the first four errors listed in Table II are only accurate 
for lensing results in a slightly narrower range: < z < 2, 
0.6 < er 8 < 1.0 and 0.1 < fi m0 < 0.5. 



Type 



Error in /i 



l/r, 



Weak-lensing approximation 
Halo model 
Mass function 
Concentration parameters 
Our fitting functions 



5% 
5% 
3% 
3% 
4% 



Total error (without fitting error) 
Grand total error 



~87o 

±10% 



Table II. Summary of the errors for /u n intrinsic in the ap- 
proach of this paper, for n = {2, 3}. For n — 4, the first four 
errors are larger, giving a total of the order of 10 — 15% 
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